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Abstract 

We discuss cross sections for hadron + hadron — > 2 jets + anything in 
which three jet variables are measured. Such cross sections are useful espe- 
cially for determining parton distributions. We define a new cross section 
da/dXA dXs drj^ for which the perturbation theory is nicely behaved even in 
the kinematic regime where the parton distributions are probed at large mo- 
mentum fractions. The cross section da/dETi drji dr]2, which has been used 
in the past, is not so well behaved in this region. We calculate these cross 
sections at order in QCD. 
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Typeset using REVT^ 



The experimental investigation of jet production in high energy hadron colhsions pro- 
vides a direct view of the underlying process, parton-parton scattering. It thus provides an 
opportunity to test quantum chromodynamics (QCD) in some detail. The one jet inclusive 
cross section da/dEx, where is the transverse energy of the jet, provides an excellent 
probe of any possible breakdown of QCD at short distances The two jet inclusive cross 
section da/dMjj, where Mjj is the invariant mass of the two jet system, plays a similar role, 
and also tests for possible resonances that might be produced in parton-parton scattering 
0. The two jet angular distribution [da/dMjj drj^]/[da / dMjj], where rj^ is half the rapidity 
difference between the two jets, probes the angular dependence of the Feynman diagrams 
for parton-parton scattering [Q]. These tests have confirmed QCD quite convincingly in a 
transverse momentum range from 30 GeV all the way to 400 GeV. 

It is possible to study inclusive two jet production in even more detail. In a Born level 
description, the physical process is the scattering of two partons. Each of the two outgoing 
partons is described by three variables. However, transverse momentum conservation elimi- 
nates two of the six total variables, while symmetry about the beam axis makes one variable 
superfluous. Thus there are three independent variables. So at the Born level one can define, 
at most, a triply differential cross section for two jet production. For instance, the rapidities 
rji and 772 of the two partons and the transverse energy E^i of one of them can serve as the 
three independent variables, leading to a cross section da / dETidrjidrj2. This choice of the 
three variables, however, is not unique and one may consider other sets of three variables to 
describe the two outgoing partons and thus other triply differential cross sections. One must 
also specify how the cross section definition is extended from the case of two partons to the 
real case of many particles. In this paper we discuss triply differential two-jet cross sections 
with the goal of using such a cross section to help determine parton distribution functions, 
particularly at large momentum fractions. We attempt to define the cross section such that 
the next-to-leading-order contributions are small compared to the Born-level contributions 
over the entire allowed phase space. 

We begin with the cross section da/dExi dr]i dri2, which has been measured by both the 
CDF and DO |P groups at Fermilab and calculated at order by Giele, Glover, and 
Kosower 0. We define jets according to the standard cone algorithm |^, supplemented 
by certain algorithms for dealing with overlapping jet cones Each jet is labeled with 

variables Et, rj and 0. Here Et is the sum of the absolute values of the transverse momenta 
of all the particles in the jet cone. (CDF adopts a slightly different definition of Et 
The variables 77 and are the ii^T-weighted averages of the rapidities and azimuthal angles 
of all the particles or calorimeter towers in the jet cone. In each event, we pick the two jets 
with the highest transverse energies. One of the two jets is defined to be the trigger jet, 
with transverse energy and rapidity Eti, rji. Then 172 is the rapidity of the other jet. Since 
there are two ways to choose which jet is the trigger jet, each event contributes to two bins 
of the cross section. 

The agreement between theory and experiment for da /dExi dr]i dri2 is satisfactory, given 
the experimental and theoretical errors However, as we shall see, the cross section 

has large theoretical errors in certain regions of the jet-variable space. In Fig. we show 
the cross section at rji = and rj2 = 2 plotted against Eti- We plot both the Born cross 
section and the full order cross section as bands whose widths refiect the theoretical er- 
ror. The calculated cross section depends on an MS renormalization scale, iiuv-, describing 
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FIG. 1. The cross section da/dETi drji dr]2 as defined in the text plotted against Eti for 771 = 
and rj2 = 2. Dark band: full order cross section. Light band: Born cross section (order a^). The 
bands reflect the variation of the calculated cross section with the renormalization and factorization 
scales. We use the CTEQ2 parton distributions 

the regulation of ultraviolet divergences, and an MS factorization scale, fxco, describing the 
regulation of coUinear divergences. We define these parameters to be fiuv = Auv Eti/^, 
lico = We choose {Auv , ^co) = (1, 1) as the standard value. Then we deter- 

mine the bands by calculating the cross sections for (log(A[/y), log(y4co)) chosen at eight 
points around the circle \o^{Auv) + \o^{Aco) = 21og^(2) = 0.96. The upper and lower 
edges of the error bands correspond to the maximum and minimum, respectively, of the 
cross section calculated in this range of (A^/y, Acq) values. Such an error band for the order 

cross section constitutes an estimate of the error induced by omitting order a^"*"^ and 
higher contributions to the cross section. 

The results shown in Fig. |l| indicate that the theory is well behaved in the central 
region, 50 GeV < Eti < 140 GeV. However, the theory is not well behaved in the region 
140 GeV < Eti- The corrections to the Born cross section are large and so are the 
estimated errors. The reason is simple p . The allowed kinematic region for the Born cross 
section for rji = and 172 = 2 is Eti < 215 GeV. Beyond this value of Eti, the momentum 
fraction of the incoming quark from hadron A becomes larger than 1. Although the Born 
cross section is zero for Eti > 215 GeV, the full cross section is not, as we will discuss 
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below. Thus for Eti > 215 GeV, and even for some interval of Eti below this value, the 
order cross section is effectively a leading order prediction. Generally speaking, leading 
order cross sections have large estimated errors. This case is no exception. The estimated 
error is greater than ±15% for Ej^i > 140 GeV. At the edge of our plot, the estimated error 
has reached ±35% and is still growing with Eti- 

The following general, if somewhat abstract, formalism will allow us to understand how 
the allowed region at order is related to the definition of the jet cross section. In an 
order calculation, one can have up to N partons in the final state. These partons can 
be described by 3A^ — 3 parameters Vt in a space Vat. The definition of a specific triply 
differential cross section, i.e. of a specific set of three jet-variables, may be understood as 
a set of maps of the A^-parton space Vat into a three dimensional space Vo of measured 
jet- variables. Equivalently, we can view the definition as mapping the Vn into V2, followed 
by a map of V2 into Vq. Thus a point v in Vat is mapped into a point Mn^v) in V2. Now 
there are many possible sets of maps A4 corresponding to different choices of the three 
variables. The only restriction comes from the requirement of infrared safety. When, in an 
A^-parton configuration, one parton becomes soft or coUinear to one of the beams or two 
partons become collinear, then that configuration must map to the same point in V2 as the 
physically equivalent (A^ — l)-parton configuration. This will ensure the cancellation of the 
relevant singularities in the perturbation theory. 

Consider, now, the physically allowed region in Vat. This is the region determined 
by Xyi < 1 and < 1, where the incoming parton momentum fractions are 

XA = j:^e^"^ XB = Y.^^-'"- (1) 

This region is mapped into a region M.]^[An) in V2 . How does M.n{,-An) compare to ^2? 
The infrared safety condition implies that 

Mn{An) ^ A2. (2) 

To see this, let t; be a two parton point in the allowed region A2. Then consider the 
physically equivalent point v in Vn obtained by adding N — 2 zero momentum partons to 
V. The infrared safety of the map implies that J^^iy) = v. Since v G ^a^, Eq- O) follows. 

While Eq. (^) holds for any infrared safe definition of a triply differential cross section, 
the answer to the question of whether J^n{An) is equal to A2 or bigger than A2 depends on 
which definition one chooses. Consider the particular cross section da/dE^i drji dri2 defined 
earher, and denote the corresponding maps as J^\. It is easy to see that the region Jv[^j^[An) 
is bigger than A2- For instance, the three parton final state with 

Eti = 220 GeV, rji = 0, 0i = 0, 

Et2 = 120 GeV, r]2 = 2, 02 = vr, (3) 
Et3 = 100 GeV, ?73 = 0, 03 = TT 

has xa = 0.67, xb = 0.19. This kinematically allowed state is mapped to the two parton 
final state with 
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Eti = 220 GeV, r/i = 0, 0i = 0, 
Et2 = 220 GeV, r/2 = 2, 02 = vr. 



(4) 



This state has xb = 0.14 but = 1.03, so it is not in the kinematically allowed region A2 
and M]yiAN) D A- 

The results depicted in Fig. Q suggest that the cross section da / dExi drji drj2 provides a 
useful tool for exploring QCD and testing parton distributions, but that there are difficulties 
resulting from the fact that the allowed regions J^n{An) at order are larger than A2, the 
region allowed at the Born level. Now the question arises, can we define a triply differential 
jet cross section such that the maps of all of the allowed regions AiN^AN) equal ^2? In 
this case the perturbative result can be well behaved over the entire allowed region. The 
answer is that there are many such solutions and one particular choice stands out as being 
particularly simple. 

We define a cross section da / dXAdXs drj^, as follows. Define jets according to the stan- 
dard cone algorithm [[^. In each event, pick the two jets with the highest transverse energies 
Et. Let 

^* = 1^1 - ^2|/2, (5) 

where 771 and 772 are the rapidities of these two leading jets. (That is, i] is the i?r-weighted 
average of the rapidities of all the particles or calorimeter towers in the jet cone 0.) Let 

= E^e-, = E^e-. (6) 

iejets V iGjets v'^ 

The sum here runs over all of the particles (or calorimeter towers) in either of the jet cones. 

Let us call the maps corresponding to this definition M^}^. It is easy to see that the 
allowed regions (^at) are all the same. For an allowed A^-parton configuration, the true 
momentum fractions Xa and x^, as given in Eq. (|1]), are less than 1. From Eq. (|^) we obtain 

Xa < Xa, Xb < xb- (7) 

Thus an allowed A^-parton state must have Xa < 1 and Xb < 1. In the two parton final 
state with the same jet variables, the momentum fractions of the final state partons are 
precisely Xa and Xb- Thus an allowed A^-parton final state is mapped into an allowed 
two-parton final state. That is, Ai^^{AN) ^ ^2- Since we already know from Eq. (^) that 
M^^{An) 5 A2, we have 

M'^.iAN) = A2. (8) 

The numerical example used earlier illustrates the point. We considered the allowed 
three parton configuration with parameters v given by Eq. (|^). With the new definition, 
this state has jet variables Xa = 0.61, Xb = 0.13, r/* = 1. The corresponding two parton 
state V = A4f{v) is 

Eti = 166 GeV, r/i = -0.23, 0i = 0, 

Et2 = 166 GeV, r?2 = 1.77, 02 = vr (9) 
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with xa = Xa = 0.61 and xb = Xb = 0.13. Thus with the new definition the two parton 
state is kinematically allowed. 

The variables Xa and Xb are conceptually straightforward. To calculate Xa-, for in- 
stance, we sum of the plus-components kf = {k^ + kf)/\/2 of the four-momenta of the 
particles inside both jet cones and divide by the plus-component of hadron A's momentum. 
On the experimental side, we note that, for large rapidity rji, ETiQ^pirji) is approximately 
twice the energy of particle i. Since calorimeters directly measure energy, the spreading 
of a large rapidity jet over a region of rapidity in the detector should not much affect the 
measurement of Xa- 

The cross section da / dXAdXB drj^ takes a very simple form at the Born level: 

X! XAfa/A{.XA, I^a) XBfb/BiXs, fJ'B) Ol^^fluv) Hab{r]*)- 
a,b 

Here fa/A{XA, f^A) is the parton distribution function for finding a parton of type a in hadron 
A, evaluated at some scale fiA', fb/BiXs, ^j^b) is the corresponding function for hadron B. 
There are two powers of as evaluated at a scale fiuv- Finally there is a function Habij]^) that 
gives the angular distribution for 2 parton 2 parton scattering at cm. energy squared s: 

da{a + h —> a' + h') ^1 . . 

ISi. = T 

The incoming partons are of type a, h here and we sum over types of outgoing partons. 
We see that the variables Xa, Xb, and rj^ are well adapted to probing the factors in the 
Born-level cross section, especially the parton distribution functions. 

What of the scales /i? Some sensible choices are suggested by the kinematics of parton- 
parton scattering at the Born level, as seen in the parton-parton cm. frame. One choice 
might be half the energy of each of the outgoing partons, E/2 = [XaXbs]^^"^ /A, while 
another might be half the transverse momentum of each of the outgoing partons, Pt/2 = 
[X^Xbs]^/^/[4 cosh(?7*)]. In our next-to-leading order calculation, we use a compromise 
between these two approaches: 



fiuv = Ajjv \/X^XBs/[4cosh(0.7?7j], 



/iA = Aco V^A^ijs/[4cosh([l - XA]r].)l 



fiB = Aco ^JXAXBs/[Acosh{[l - Xb]v*)]- (12) 

Here the factors Auv and Aco are adjustable. Any choices of these ^'s that are of order 1 
would be reasonable. 

In Fig. 1^, we show the cross section da/dXAdXB drj^ at Xa = 0.1 and rj^ = plotted 
against Xb- We plot both the Born cross section and the full order cross section, in 
each band whose width refiects the theoretical error as determined by choosing 

(log(A[/y), log(y4co)) at eight points around the circle \og^{Auv) + log^(^co) = 21og^(2) = 
0.96, just as in Fig. |I|. 
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FIG. 2. The cross section da/dXAdXB dr]:^ as defined in the text plotted against Xb for 
Xa = 0.1 and r/^, = 0. Dark band: full order cross section. Light band: Born cross section 
(order a^). The bands reflect the variation of the calculated cross section with the renormalization 
and factorization scales. We use the CTEQ2 parton distributions Q. 



The results shown in Fig. ^ indicate that the theory is well behaved in the entire Xb 
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range shown, even out to quite large Xb- The corrections are not large and the a 
calculation has small estimated errors from higher order contributions, about 10%. 

Our calculations in this paper are based on the subtraction algorithm and corresponding 
computer code for calculating next-to-leading order jet cross sections described in Ref. JlOl . 



The desired definitions of jets and jet variables are inserted into certain small subroutines 
in the code. A new feature compared to previous versions of the code is that the program 
directly calculates the so-called i^'-factor, the ratio of the full order cross section to the 
Born cross section. This has technical advantages stemming from the fact that the K- 
factor is normally nearly constant as a function of the three jet variables, at least for an 
appropriately defined cross section as here. 

We conclude with a caveat. We have seen that the cross section da/dX^dXB drj^, is 
reasonably well behaved for small 1 — Xa or 1 — Xb- Nevertheless, the example of the Drell- 
Yan cross section teaches us that the perturbation series will contain terms with factors of 
log(l — Xa) and log(l — Xb)- A summation of these logarithmic terms will be useful, and 
perhaps necessary, in order to explore with precision the small 1 — Xa and 1 — Xb regions. 
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Substantial progress along these lines has been made in the Drell-Yan case [ITT 
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